arXiv:1505.01685vl [cond-mat.str-el] 7 May 2015 


Flat-band ferromagnetism in a topological Hubbard model 

R. L. Doretto 

Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, 13083-859 Campinas, SP, Brazil 

M. O. Goerbig 

Laboratoire de Physique des Solides, CNRS UMR 8502, 

Universite Paris-Sud, F-91405 Orsay Cedex, France 
(Dated: May 8, 2015) 

We study the flat-band ferromagnetic phase of a topological Hubbard model within a bosonization 
formalism and, in particular, determine the spin-wave excitation spectrum. We consider a square 
lattice Hubbard model at 1/4-filling whose free-electron term is the yr-flux model with topologically 
nontrivial and nearly flat energy bands. The electron spin is introduced such that the model either 
explicitly breaks time-reversal symmetry (correlated flat-band Chern insulator) or is invariant under 
time-reversal symmetry (correlated flat-band Z 2 topological insulator). We generalize for flat-band 
Chern and topological insulators the bosonization formalism [Phys. Rev. B 71, 045339 (2005)] 
previously developed for the two-dimensional electron gas in a uniform and perpendicular magnetic 
field at Ailing factor n = 1. We show that, within the bosonization scheme, the topological Hubbard 
model is mapped into an effective interacting boson model. We consider the boson model at the 
harmonic approximation and show that, for the correlated Chern insulator, the spin-wave excitation 
spectrum is gapless while, for the correlated topological insulator, gapped. We briefly comment on 
the possible effects of the boson-boson (spin-wave-spin-wave) coupling. 

PACS numbers: 71.10.Fd, 73.43.Cd, 73.43.Lp 


I. INTRODUCTION 

Electronic bands with non-zero Chern numbers are at 
the origin of a large variety of topological phenomena 
in condensed-matter systemsIn a pioneering work in 
1988, Haldane showed that a two-dimensional graphene¬ 
like lattice model with broken time-reversal symmetry 
can exhibit an integer quantum Hall effect (IQHE) with¬ 
out an external magnetic field.l^ Later, in 2005, Kane 
and Mele generalized this model to restore time-reversal 
symmetry with the help of the natural spi n de gree of free¬ 
dom in graphene with spin-orbit coupling.®^ The lowest- 
energy spin bands in this model carry non-zero but op¬ 
posite Chern numbers that result in the quantum spin 
Hall effect (QSHE), which manifests itself in a quan¬ 
tized conductance associated with the transverse spin 
current. Whereas the intrinsic spin-orbit coupling is too 
small in graphene to reveal the effect, the QSHE was 
later predicted^ to occur and measurecP in HgTe/CdTe 
quantum wells. 

The presence of an IQHE in a band with a non-zero 
Chern number indicates a certain similarity between the 
band and the Landau level of the two-dimensional elec¬ 
tron gas (2DEG) in a strong magnetic held. However, 
in contrast to the latter the energy bands obtained in 
tight-binding models have usually a non-negligible dis¬ 
persion. In order to investigate in further details the 
relation between Landau levels and energy bands with 
non-zero Chern numbers, special effort has recently been 
invested into the engineering of flat bands in specially de¬ 
signed tight-binding models.^^®^ If these bands are par¬ 
tially fllled and if electron-electron interactions are taken 
into account, one would then expect correlation effects 


similar to the fractional quantum Hall effect (EQHE). 
The effect, also called fractional Chern i nsulat or, was 
later corroborated within numerical studie^^^ (for re¬ 
cent reviews, see ISHIil). 

The analogy between flat bands with non-zero Chern 
number and Landau levels can be pushed further when 
the internal spin degree of freedom is taken into ac¬ 
count. Indeed, when there are as many electrons as flux 
quanta threading the 2DEG, the spins are spontaneously 
aligned and form a ferromagnetic state (quantum Hall 
ferromagnet) in order to minimize the electron-electron 
repulsion.E^ This situation corresponds to half-filling in 
a lattice model if only the two lowest (spin) bands are 
taken into account. However, in contrast to Landau lev¬ 
els, where time-reversal symmetry is broken by the exter¬ 
nal magnetic field and the Landau levels occur merely in 
two spin copies, the situation is more involved in lattice 
models. One needs to distinguish two generic situations. 
In the first one, time-reversal symmetry is preserved such 
that spin and orbital degrees of freedom are coupled. In 
this case, possible spin excitations are described in the 
framework of rather unusual commutation relations for 
the spin-density operator and the resulting ferromag¬ 
netic state is expected to respect the underlying Z 2 sym¬ 
metry of topological insulators.® This type of ferromag¬ 
netism has recently been investigated within numerical 
exact-diagonalization studies by Neupert et al, who find 
a gapped Ising ferromagnetic ground state.l^^The second 
situation arises when time-reversal symmetry is broken 
on the level of the lattice model, and where the Chern 
bands occur in two spin copies without any spin-orbit 
structure. 

In the present paper, we investigate ferromagnetism in 
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both situations, within a specially adapted tight-binding 
model on a square lattice with on-site Hubbard repul¬ 
sion that is a generalization of the model originally pre¬ 
sented in Ref. uni The tight-binding model, in the ab¬ 
sence of interactions, bears a staggered tt-Aux phase for 
each spin component, and time-reversal symmetry de¬ 
termines whether the two spin species experience either 
the same or an opposite flux per plaquette. If time- 
reversal symmetry is broken, one is confronted with a 
correlated flat-band Chern insulator, whereas one finds a 
correlated flat-band Z 2 topological insulator in the case of 
preserved time-reversal symmetry. For both situations, 
we investigate the ferromagnetic state at quarter-filling 
of the lattice that corresponds to half-filling of the two 
lowest energy bands. In order to investigate its stabil¬ 
ity and collective spin-wave excitations, we construct a 
nonpertubartive bosonization scheme similar to the one 
proposed in Ref. [19] to describe the 2DEG at Ailing fac¬ 
tor V = 1. Such a formalism was also applied to study 
quantum Hall ferromagnetic phases realized in graphene 
at filling factors 1 / = 0 and ly = ± Poland to describe the 
Bose-Einstein condensate of magnetic excitons realized 
in a bilayer quantum Hall system at total filling factor 

= imm 


A. Overview of the results 

We show that the bosonization scheme,^ originally de¬ 
veloped for the 2DEG at Ailing factor 1 ^ = 1, can be gen¬ 
eralized for lattice models that describe flat-band Chern 
and Z 2 topological insulators. 

For both correlated Chern and Z 2 topological insulators 
described above, we map the interacting fermion model 
to an effective interacting boson model. We consider the 
effective boson model in the harmonic approximation and 
show that the ground state is indeed given by a spin po¬ 
larized (ferromagnetic) state. Our main results are in 
fact the analytical calculation of the spin-wave excita¬ 
tion spectra of both ferromagnets (Figs. |^and[^. The 
spin-wave dispersion relation corresponds to the energy 
of the bosons at the harmonic approximation. Due to 
the bipartite nature of the underlying square lattice, we 
identify two types of collective spin-wave excitations. We 
And that the correlated flat-band Chern insulator has one 
gapless and one gapped spin-wave excitation branches 
(Fig. 1^ while the correlated flat-band topological insu¬ 
lator, two gapped ones (Fig. [^. For the latter, the exci¬ 
tation gap we obtain at zero-wave vector coincides with 
the result numerically calculated by Neupert et 


B. Outline 

Our paper is organized as follows. Section [H] intro¬ 
duces the basic tight-binding model (the spinfull square 
lattice TT-flux model), which is discussed in view of the 
role played by time-reversal symmetry. We discuss the 


flat-band ferromagnetic phase of a correlated Chern in¬ 
sulator with broken time-reversal symmetry in Sec. jllH 
whereas the more involved case of a model with underly¬ 
ing time-reversal symmetry (correlated topological insu¬ 
lator) is presented in Sec. IV For both cases, the particu¬ 
larities of the associated lattice models are first discussed 
in Secs. [HI Aj and jlV Aj before we present the details of 
the bosonization schemes in Sec. jlUBj and Sec. jlVB] re¬ 
spectively. The bosonization formalism is then applied in 
Secs. [HIC] and |lVC| to study the flat-band ferromagnetic 
phases obtained in the presence of a strong on-site Hub¬ 
bard repulsion term in the respective models. We com¬ 
ment on possible extensions of the bosonization scheme 
and the effects of the boson-boson (spin-wave-spin-wave) 
coupling in Sec. |V] and, in Sec. lYg we provide a brief 
summary of our findings. Technical details of the two 
bosonization schemes, as well as a more detailed analysis 
of time-reversal symmetry, are delegated to the (three) 
Appendices. 


II. TIGHT-BINDING MODELS WITH FLAT 
TOPOLOGICAL BANDS 


Before a detailed analysis of the different ferromag¬ 
netic states in flat-band Chern insulators with broken 
time-reversal symmetry and Z 2 topological insulators, we 
present here a common tight-binding model that provides 
the different flat bands. 

Let us consider = N free spin-1/2 electrons hopping 
on a bipartite square lattice where both sublattices A and 
B have each = Ng = N sites. The Hamiltonian of 
the system is given by the tight-binding model 

iG.A, n, (T 

+ 'y ^ ^Ci-i-5 Q o- + H.C.^ . (1) 

2 , 6, a, (T 

Here (ciao-) creates (destroys) a spin a =t,4- 

electron on site i of the sublattice a = A,B. The 
spin-dependent nearest-neighbor j-,cr and next-nearest- 
neighbor Xi j hopping energies are respectively given by 

[Fig-lIJa)] ’ 

ti,i+n,a = h exp [i(-l)”7(cr)7r/4], i G A, (2) 


_ (-{-iyt2, iGA, 

iGB, 


( 3 ) 


where ti and ^2 are positive real quantities. The indices 
n = 1,2,3,4 correspond to the nearest-neighbor vectors 
[Fig.[3b)] 

Ti = -Ta = - (ai - a 2 ) = - y), 

( 4 ) 

1 / X 0 -,- 

T2 = -T4 = - (ai -g a 2 ) = 2 ^^ + v)^ 
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FIG. 1. (Color online) (a) Schematic representation of the hopping term 0. The nearest-neighbor hopping energies (solid 
black lines) are equal to t\ exp (z7r/4) (spin up electrons) in the direction of the arrows and the next-nearest-neighbor hopping 
energies are equal to -\-t 2 (dashed red lines) and —t 2 (thin blue lines). Black and open circles indicate sites of the A and B 
sublattices, respectively, (b) The nearest-neighbor vectors Ti [read arrows, Eq. @] and the primitive vectors ai = ax and 
bl 2 = ay [green arrows, Eq. 0]. The next-nearest neighbor distance a is set to one. (c) Brillouin zone. Here X = (tt, 0), 
Y = (0,7r), and M = ('7r,7r). 


and S is either 1 or 2 corresponding to the next-nearest- 
neighbor vectors 

= ai = ax, S 2 = a .2 = ay. (5) 

In the remainder of this paper, we set the next-nearest- 
neighbor distance a = 1. The spin-dependent phases 
7 (cr) reflect that time-reversal symmetry is preserved for 
7 (t) = — 7 ( 4 .) = 1, whereas it is broken if 7 (t) = 7 ( 4 .) = I 
(see Appendix for details). The tight-binding model 
([^ is a generalization for the case of electrons with spin of 
the TT-flux model discussed in Ref. uni Since the nearest- 
neighbor hopping energy is complex, each electron 
acquires a phase tt as it hops around a plaquette in the 
direction of the arrows indicated in Fig. 3a)- There¬ 
fore, Hq describes noninteracting electrons hopping in a 
square lattice in the presence of a fictitious staggered ±7r 
flux pattern.l^ For the time-reversal-symmetric model 
[ 7 (t) = — 7 ( 4 -) = 1 ] the flux experienced by electrons 
of opposite spin is opposite, whereas it is the same for 
7 (t) = 7 ( 4 -) = he., in the case of broken time-reversal 
symmetry. 

After Fourier transformation, i.e., introducing 

exp(-fk • ( 6 ) 

kGBZ 


is a four-component spinor. Furthermore, 

hi, = So.kTo + Bk • T (10) 

is a 2 X 2 matrix where tq =12x2 is the identity matrix 
and f = (ri, T 2 , T 3 ) is a vector whose components are the 
Pauli matrices. Finally, i?o.k and the components of the 
vector Bk = (Ri,k, R 2 ,k, Ba.k) are given by the functions 

Ro,k = 0, 

iJi k = 2V2ti cos ^ cos 

R 2 ,k = 2V2ti^{a) sin y sin 

Ra.k = 2^2 (cos kx — cos ky) . ( 11 ) 

Again, the factor 7 ( 0 -) indicates whether time-reversal 
symmetry is broken or not. Whereas the time-reversal- 
symmetric model [ 7 (t) = — 7 ( 4 -) = 1 ] has the usual prop¬ 
erty hi = (hi^)*, we find that the two components of 
the Hamiltonian (|^ are identical for all wave vectors, 
hi = hi, in the case of 7 (t) = 7 (i) = 1 . 

A. Symmetries of the 7r-flux model: spin rotation 


with the momentum sum running over the BZ associated 
with the underlying Bravais lattice [see Fig. 3 c)], it is 
possible to show that the hopping term 0 assumes the 
form 

^ 0 = Y. 'f'k^k'I'k, (7) 

kGBZ 


where 

nk = 

is a 4 X 4 matrix and 


'Sl = 


0 hi 


'^kAt '^kSt '^kAi CkBi 


( 8 ) 

(9) 


The above discussion and the role of time-reversal sym¬ 
metry allow us to investigate certain properties of the 
supposed flat-band ferromagnetic states from a pure sym¬ 
metry point of view. In this section, we discuss the be¬ 
havior of the noninteracting fermion model (3 under spin 
rotation. Some further considerations about the behavior 
of 3 under time-reversal are present in Appendix [3 
The Hamiltonian 3 can indeed be written as 

( 12 ) 

k,(7 a,b 


Comparing Eq. © with Eqs. 3 and ( [To| , we see that, 
for the Chern insulator with broken time-reversal sym- 
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metry {hl = h^), 

E^^{k) = El^{k) and Ef{k) = Ef{k) (13) 

with o 7 ^ 6 , while for the topological insulator with time- 
reversal symmetry {h^ = we have 

E^^{k) = El^{k) and (k) 7 ^ (k) (14) 

with a ^ b. 

The total spin operator reads 

s = Sj(j — — Cp^^CTo-.o-'Cpacr' (15) 

i,a a,p 


where is the spin operator at site i of the sublattice 
a and a = (cti, ( 72 , us) is a vector of Pauli matrices as¬ 
sociated with the physical spin [see Eq. (25) below]. It 
is possible to show that the following commutation rela¬ 
tions hold: 


[no,s^] = o 


[^0, = E E [^r(k) - Eri^)] 4, 

k a 

+ [Ef^ik)-Ef^{k)]cl^^c^^Bi 

+ [Ef^(k)-Ef^(k)]cEtCk.A4 


L 4- 5 


[Ho,S-] = [S+,Hor. 


(16) 


Since the Chern insulator is characterized by the 
Hamiltonian ([^ with two identical components for the 
two spin orientations [see Eq. (13)], one immediately re¬ 
alizes that all components of the total spin operator (15) 
commute with the Hamiltonian ([^, i.e., the Hamiltonian 
has SU(2) spin rotation symmetry. In the case of a ferro¬ 
magnetic ground state that breaks spin rotation symme¬ 
try, all equivalent states can thus be obtained by global 
rotations generated by the total spin operator (15), and 


one would therefore expect the presence of a Goldstone 
mode in the spin-wave excitation spectrum. We show 
explicitly in the following section the existence of such a 
mode. 

Concerning a topological insulator described by the 
Hamiltonian <§ with coefficients obeying Eq. ( [l4| , one 
can easily see that the SU(2) spin rotation symmetry is 
now explicitly broken to U(l), i.e., the Hamiltonian is 
invariant under spin rotations around the z-axis. There¬ 
fore, the presence of a Goldstone mode depends on the 
ground-state spin polarization: whereas one would ex¬ 
pect a superfluid-type mode for an easy-plane ferromag¬ 
netic state, where the polarization is oriented in the xy- 
plane, an easy-axis ferromagnet with a polarization along 
the z-direction would not display a Goldstone mode since 
the ground state preserves the U(l) spin rotation sym¬ 
metry of the Hamiltonian. In the following, we show that 
the latter is the case for a topological insulator and that 
all collective excitations are indeed gapped. 


III. FERROMAGNETISM IN A FLAT-BAND 
CHERN INSULATOR 

In this section, we discuss the flat-band ferromagnetic 
phase of a topological Hubbard model that explicitly 
breaks time-reversal symmetry. We start by discussing 
the free electron term of the model. 


A. Square lattice tt-Hux model with broken 
time-reversal symmetry 


The noninteracting Hamiltonian (0 with 7 (t) = 7(4- 
) = 1 can be diagonalized with the aid of the Bogoliubov 
transformation 


'"kScr “ ^k(^ko. ~ Wk'-kcr, (17) 


where the Bogoliubov coefhcients Uk and Uk are given 
by Eq. (Bl). After diagonalization, the Hamiltonian 0 
assumes the form [see Eqs. (|B1[)-(B3) for details] 


^ 0 = E ^Ic^'k^k, (18) 


where the 4x4 matrix reads 


n^i.= 



with 


h 


/ _ 

k — 


‘^(i,k 0 \ 

0 “^Cjk J 


(19) 


( 20 ) 


being a 2 x 2 diagonal matrix whose elements are the 
upper band d (- 1 - sign) and the lower one c (— sign), 


Wd/c.k = -So.k ± |Bk|, (21) 

and the new four component spinor <i>]{^ is given by 

The band structure of Hq is made out of four bands: 


dt, di, ct, ci, 

with the c and d bands doubly degenerate in the spin 
degree of freedom as expected from the discussion of the 
previous section. Figure shows the energy of the c 
and d bands [Eq. @] for different values of the ratio 
t 2 /ti. Note that the spectrum is gapless for t 2 = 0 
and, as ^2 increases, it acquires a finite energy gap 
A = min(a;d,k) — max(a;c,k)- Indeed, as shown in Fig. 

A linearly increases with t 2 for t 2 < 0.5 H and then it 
saturates at A = 4ti for larger values of ^ 2 - The widths 
HA/d = max(a;c/d_k) ~ niiii(wc/d,k) of the c and d bands 
also change with ^ 2 - In particular, as illustrated in Fig.[^ 
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FIG. 2. (Color online) Band structure of the noninteract¬ 
ing hopping term T-Lq [Eq. ( |21[ ) in units of ti] along paths in 
the Brillouin zone [Fig. [^c)J for different values of the next- 
nearest-neighbor hopping energy amplitude t 2 . Solid green 
and dashed blue lines correspond respectively to the c and d 
bands. 


the flatness ratio fc = A/Wc of the c-bancP^ is maxi¬ 
mum (fc = 4.83) for 0.5 < 12/^1 < l/v^- Such a range 
includes the configuration <2 = ii/'s/S discussed by Neu- 
pert et 

An interesting aspect of the tight-binding model Hq 
is that its energy bands are topologically nontrivial. In¬ 
deed, it is possible to show that the Chern numbers of the 
c and d bands can be written in terms of t he coeffic ients 
Bi,k [Eq- ( [lT| )] and that they are finite, i.e., b l ^O0 | 2 ll 

= [ d^k Bk ■ {dkjk X dk^Bk) = Tl, (23) 

47r Jbz 


with Bk = Bk/|Bk|, regardless the value of ^2 > 0. 
Therefore, the square lattice 7r-flux model Q with bro¬ 
ken time-reversal symmetry is an example of a Chern 
insulator. As discussed in the Introduction, the system 
should exhibit an IQHE when a certain number of energy 
bands are completely filled. 

In the following, we focus on the nearly flat-band limit 
of 'Ho-, in particular, we consider ^2 = ti/'s/S. Note that 
since each of the c a and d a bands have = Ng = N 
available states and = N, the d band is empty while 
the c one is half-filled, i.e., we have 1/4-filling includ¬ 
ing all four bands. In the next section, we introduce a 
bosonization scheme to describe such a flat-band Chern 
insulator. 



FIG. 3. (Golor online) Flatness ratio A/W of the c-band 
and band gap A as a function of the next-nearest-neighbor 
hopping energy amplitude t 2 for the noninteracting hopping 
term Ho- 


B. Bosonization scheme for a fiat-band Chern 
insulator 

In order to study the flat-band ferromagnetic phase of 
a correlated Chern insulator, we introduce a bosonization 
scheme similar to that proposed in Ref. [Tni for the 2DEG 
at zz = 1. Following the lines of the bosonization method 
for one-dimensional fermion systems,!^ the idea of the 
formalisirfi^ is to define boson operators in terms of the 
lowest energy neutral excitations of the original fermionic 
system and then map the interacting electronic model to 
an effective bosonic one. In the following, we outline 
the bosonization scheme for flat-band Chern insulators. 
The (two) approximations involved and the differences 
between the scheme on a lattice and the one developed 
for the 2DEG at 1 / = 1 are discussed. Further details can 
be found in Appendix [B} 

Let us consider the tight-binding model 0 at 1/4- 
filling. As mentioned in the previous section, in this case 
the highest energy d bands are completely empty while 
the lowest energy c ones are half-filled, see Fig. In 
particular, let us assume that the c f band is completely 
filled while the c one is completely empty, i.e., we con¬ 
sider that the ground state of the free electron model Q 
is completely spin polarized, 

|FM) = n c|,^|0). (24) 

kGBZ 

The state |FM) is the reference state of the bosoniza¬ 
tion method and belongs to the ground-state manifold 
of ferromagnetic states that need to be considered when 
electronic interactions are taken into account. Since the 
lowest energy c bands are separated from the highest 
energy d ones by an energy gap and the former bands 
are partially filled, the lowest-energy neutral excitations 
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are given by particle-hole pairs (spin flips) within the 
c bands. Therefore, in the following, we neglect the d 
bands, i.e., we restrict the Hilbert space to the subspace 
spanned by the lowest energy c bands. 

The spin operator at site i of the sublattice a is defined 
as 


where the gQ;(p,q) functions are defined as 


5a(p, q) = Up. (33) 


Similar considerations hold for the density operator of 
electrons with spin a at site i, which is define as 


1 


.t 


Si a — ^CJfj fj'Ci 


(25) 


in terms of the Pauli matrices defined above. The Fourier 
transform of the components of Si a read 


Pi a G a G^i o. G ■ 

Its Fourier transform is given by 


^ X! exp(zq • R*)/3a„(q) (35) 


(34) 


with X = x,y, z. In particular, using Eq. 
to show that 


SL = exp(iq • R*)5'q,a (26) 

(|^, it is possible 

^q,a ^ g g t*"P ° 

P 

^ ( ^p—g a.l,*"P ° T’ (27) 

P 

where ± . With the help of the Bo- 

goliubov transformation ( |I7[ ), one can express the spin 
operators in terms of the c and d fermion operators. 
For instance, 

(^P-q^p'^p-qt'^pf ^p-q^p''p-qt''Pf 


N 

qGBZ 

and, with the help of Eq. (§, it is easy to see that 

Pao-(q) = ^ ( Cp_q ^ jCpg o-. (36) 

p 

Finally, the projected electron density operator reads 

Pa<T(q) = Xl‘^“(P’^)4-q<TCpcT (37) 


with the Ga(p,q) functions given by Eq. (30). 

Differently from the Girvin-MacDonald-Platzman 
(GMP) algebrrP^ for electrons within the lowest Lan¬ 
dau level here, for a flat-band Ghern insulator, the al¬ 


gebra of the projected spin [Eq. (31)] and electron den¬ 
sity [Eq. (37)] operators is not closed. For instance, the 
commutator 


+ “p-q'(^pC^p-qtCpt + 'Pp-q?^pCp_q^dp;) , (28) [S^^^ , _^] = ^ (P “ q', q)5/3 (P, HOc, 


p-q-q' 


t'^pt 


where Up and Vp are the Bogoliubov coefficients (Bl|. 
Note that by projecting 5'+^ into the c bands, only the 


second term on the R.H.S. of Eq. (28) survives. Indeed, 


it is easy to see that the expressions of the projected spin 
operators read 


P 

>5q,a = 


where 


Gyi(p,q) = Up 


and Gb(p, q)=M: 


p-q“p 


(29) 


• (30) 


- 5a(p,q)5/3(p- q,q')Cp_q_q/^Cp; (38) 

cannot be expressed in terms of the projected spin and 
electron density operators [see, for instance, Eq. (37) 
from Ref. [T^ . Similar considerations hold for the com¬ 
mutators 

[Pa CT (k), , [pa o- (k) , <Sq Q,], [pa cr (k) , Pb a' (q)] > 


As discussed below, it is interesting to consider the fol¬ 
lowing linear combination of the spin operators 

(31) 

with X = x,y,z and a = 0,1. We have, for instance, 

P 


see Eqs. (B4|-(B6). We refer the reader to Refs. IT71 l?7l 
and [5S] for a discussion about the flat Berry curvature 
limit, where the GMP algebra holds for flat-band Ghern 
insulators. 

In order to define boson operators in terms of the 
fermion operators c, we consider the neutral particle-hole 
pair excitations above the ground state jFM) of the free- 
electron model (18). Such excitations, which correspond 


to spin-flips, are given by 

14/q) = JFM), 


a = 0,1, 


(39) 


where ^ is the linear combination (31) of the spin op¬ 


erators iSq ^ and S'qB- shown below, the fermionic 
representation of the spin operators 5q ^ allows us to de¬ 
fine two sets of independent boson operators. 
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The commutator (38) between the spin operators 


and S~ p differs from the usual canonical commutation 
relation between creation and annihilation boson opera¬ 
tors. However, if the number of particle-hole pair exci¬ 
tations is small, one can write 


Cp-qt^Pt « (FM|c^_q^Cpt|FM) = (5q,0, 

« (™|c^_q^Cp;|FM) = 0. (40) 


In this case, the commutator (38) acquires the form 


« dq,-q'^5a(p-q',-q')5/3(p,q')- (41) 


Moreover, it is possible to show that for the square lattice 
TT-flux model (18) the sum over momentum in the above 


equation is finite only if a = /?, see Eq. (B7) for details. 


i.e., 


consider the projected electron densi ty o perator pa^(k) 
[Eq. From Eqs. (37), ( [4^ , and (B5), one finds that 


[Pat(k).^Lq] = - 


Ga(p,k) 


p 


F, 


a,q 


5a(p k, q)Cp_k—q^Cp^. 

r -1 

Differently from the 2DEG at v = it is not possible 
to express the commutator (46) in terms of the boson 
operators 5 q, and therefore, it is not easy to determine 
the expansion of /Oa^(k) in terms of the bosons ba that 
satisfies the commutator (46). This is related to the fact 


that for the Chern insulators the algebra of the spin and 
electron density operators is not closed (see discussion 
above). We then proceed as follows: In principle, we can 
write 






(47) 


['S'q.a.'S'q,^^] « dq,_q-3“(P “ (P’^') ■ 

Therefore, as long as the number of particle-hole pair ex¬ 
citations above the reference state |FM) is small, the com¬ 
mutator (38) is approximately equal to a canonical boson 
commutation relation. In other words, in this limit, the 
lowest-energy particle-hole pair excitations can be ap¬ 
proximately treated as bosons. We then define two sets 
of independent boson operators 


9+ 

b = 

'^“■q p 

J- rv Cl 


1 




Xl5«(P’-q)4+qT^qH 


Q- 

^ Ck.q 




a,q 






(43) 


with q: = 0,1, that obey the canonical commutation re¬ 
lations 


k,q) 


[&a,k, &/3,q] = [^X’ ^kq] “ 


(44) 


Here, the gQ(p,q) functions are given by Eq. (33) and 


Interestingly, it is possible to show that the Fa,q func¬ 
tions can be explicitly written in terms of the Bi p coef¬ 
ficients, see Eq. ( |B8[ ). It is worth emphasizing that the 
boson operators ba are defined with respect to the refer¬ 
ence state |FM). 

Once the boson operators ba are defined, we can de¬ 
rive the bosonic representation of any operator O that 
is expanded in terms of the fermion operators c. As dis¬ 
cussed in Ref. mi such a procedure consists of calculating 
the commutators [O, b^, q], writing them in terms of the 
boson operators ba, and determining the action of the op¬ 
erator O in the reference state |FM). For instance, let us 


where the i7a(k, p, q) function satisfies the relation [com¬ 
pare Eqs. (@) and @1 

X 5a (p> q)-f^/3(k, p, q) = da,/34,qFa,q. (48) 


With the aid of Eq. ( |47| ) , the commutator ( |46| ) reads 

[4T(k),&l,q] = - ^XX^“(P’^)5a(p- k,q) 


k',p 


xi7/3(k',p,k-bq)6^ 


k'- 


(49) 


It is then easy to see that the expansion 

1 


Pat 


«=-xx 


_ _ F 

o 1 , ^ Q;,q 

Oi,p k' ,p,q 


Ga(p,k)5a(p - k, q) 


xi7/3(k',p,k-bq)6X-&a,c 


(50) 


of pat(k) in terms of the bosons ba satisfies the commu¬ 
tator (49). Apart from a constant (see below), Eq. (50) 
corresponds to the bosonic representation of the electron 
density operator P(j^(k). 

Although it is difficult to solve Eq. (48) and calculate 
i7a(k, p, q), it is indeed possible to determine the partic¬ 
ular value for k = q, i7Q,(q, p,q). Comparing Eqs. (45) 
and (48), we see that 

-ffa(q,p,q) = T;^5a(p-q,-q)- (5i) 

^ Q:,q 

Note that by keeping only the term k' = k 4- q in the 
momentum sum in Eq. (50) and using Eq. (51), an ap¬ 
proximated bosonic representation for the electron den¬ 
sity operator pQ^(k) follows. Similar considerations hold 
for Pa I (k). We then arrive at the following approximated 
bosonic representation for the electron density operator: 


p,,(k)« -7V4,t4.o + XX^° 

0!,/3 q 


,(k,q)6l,,, 


k-|-q“a.q- 


( 52 ) 



























Here, the first term is related to the action of patT(k) in 
the reference s tate |FM) and the Gapaa{x,y) function is 
given by Eq. (B9) in Appendix [B| Similar to the q 
function (45), the Gapacrix^y) function can also be ex¬ 
pressed in terms of the coefficients [see Eq. (BIO)]. 


In summary, we show that the bosonization formal¬ 
ism introduce in Ref. [m for the 2DEG at iz = 1 can be 
extended to flat-band Chern insulators with a half-filled 
energy band. In both cases, it is possible to show that 
the particle-hole pair excitations can be approximately 
treated as bosons as long as the number of such exci¬ 
tations is small. While for the 2DEG an exact bosonic 
representation for the electron density operator can be 
derived, here, due to the fact that the algebra of the 
spin and electron density operators is not closed, only 
an approximated bosonic representation for the electron 
density operator can be obtained. We should note that 
the approximated bosonic representation for the electron 
density operator (52) is similar to the exact one derived 


for the 2DEG at = 1 [see Eq. (27) from Ref. [19]. 


C. Topological Hubbard model I 

Let us now consider a square lattice Hubbard model 
at 1 /4-filling whose Hamiltonian is given by 


T~Lch = "Ho + Hu- 


(53) 


Here, Hq is the tight-binding model ([^ with t 2 = tily/2 
(nearly flat-band limit) and 


'^U = U'^ X/ PiapPial 
i a—A,B 


(54) 


is the one-site Hubbard term with Piap being the electron 
density operator (34) and U > 0. In momentum space, 
TLu reads 


= /5aT(-k)/3a;(k) 


N 


(55) 


with Paaiy) given by Eq. ( |36| and N being the num¬ 
ber of unit cells, as mentioned before. Since the choice 
t 2 = implies that the energy bands of "Ho have 

nonzero Ghern numbers [Eq. (|23|)], the Hamiltonian (53) 


corresponds to a topological Hubbard model. In the fol¬ 
lowing, we apply the bosonization formalism introduced 
in the previous section to study the flat-band ferromag¬ 
netic phase of the correlated Ghern insulator (53). 

We start by projecting the Hamiltonian (53) into the 
lowest energy c bands: 


Here 


^Ch —t 'H-Ch = 'Ho + 'H.IJ. 


^0 — ^ ^ ^c.q^gcr^qc 

q 


(56) 


(57) 


[see Eq. (18)] and "He/ is given by Eq. (55) with the re¬ 
placement ; 0 acr(k) —)■ pao-(k). Some issues about the rel¬ 
evant energy scales need to be emphasized: On the one 
hand, in order for the projection scheme to the lowest c 
bands to remain valid, the energy scale U must be smaller 
than the energy separation A between the bands c and 
d, see Fig. Otherwise, the on-site interaction would 
mix the different bands and it would no longer be valid 
to characterize them in terms of Ghern numbers associ¬ 
ated with the non-interacting model. On the other hand, 
we consider the on-site interaction U to be much larger 
than the bandwidth of the (almost flat) c bands, such 
that their dispersion may be neglected in the remainder 
of the section. In this sense, the c bands are reminiscent 
of the highly-degenerate flat Landau levels of a 2DEG in 
a strong magnetic field. 

Following the same procedure used to determine the 
bosonic representation of the projected electron density 
operator (37), we show that the bosonic representation 


of Ho is simply a constant Eq [recall that in the flat-band 
limit, the kinetic energy is quenched, see Eqs. (Bll)- 


(B14) for details]. The bosonic representation of the pro¬ 
jected on-site Hubbard term Hu can be easily derived by 
substituting Eq. (52) into Hu and normal ordering the 


resulting expression. We then arrive at the boson model 


Hb=Eo+H%'^ + H 


/(4) 


(58) 


where Eo = 2.44A, the quadratic boson term is given by 




a,/3 pGBZ 


and the boson-boson interaction term reads 


_(4) _ _1_ .raPa'P'A ,t h h , 

*^k.q,p a/3'.p+k"/3.q-k"a,q"a5P' 


1 


' k,q,p 


(60) 


Here 


U 


^P^ = y E^“/3a;(0,p) 


2 

U 

N 

U 


E ^^a'at(k>P)0ac'a;(k,p), (61) 


a,a' ,k 


= ^Y.^*aPati-^,C^)Ga' P'aliKp), (62) 


with the Gapaa{^,p) function given by Eq. (B9) and 
a,/3, a',/?' = 0,1. Therefore, within the bosonization 


scheme introduced in Sec. HIB the Hubbard model (53) 


is mapped into the effective interacting boson model (58). 


In order to discuss some features of the effective boson 


model (58), let us first neglected the quartic boson term 


H^g and consider only 


HB-EoHHf. (63) 

Such a lowest-order approximation is the so-called har¬ 


monic approximation. The quadratic Hamiltonian (63) 
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can be easily diagonalized, namely 


T~Lb — Eq + ^ ^ ^ ^ ( 64 ) 

/i—=h p^BZ 


where the dispersion relations of the bosons a± read 


with 


n 


±.p 


= ±Ap + 


,2 _ ,10,01 

P P P 


(65) 


e - -le™ 

P “ 2 ^ P 


Ap=2[4°-4" 


( 66 ) 


and the given by Eq. (61). 

The ground state of the Hamiltonian (64) is the 
vacuum state for the bosons a^. Since &Q,k|FM) = 
ap,k|FM) = 0, the ground state of (64) is indeed the 
spin polarized ferromagnet state (24). This result in¬ 


dicates that the topological Hubbard model (53) has a 
stable flat-band ferromagnetic phase. 

The stability of the ferromagnetic ground state is also 
corroborated by the dispersion relations fl±.p for the 
bosons a± (Fig. |^, which correspond to the spin-wave 
spectrum of the flat-band ferromagnetic ground state 
|FM). One notices that the energy scale of the excita¬ 
tions is given by the on-site repulsion energy U instead 
of the nearest-neighbor hopping energy ti since the boson 
representation of "Ho is simply a constant Eq. The exci¬ 
tation spectrum has two branches, a gapless one (H+.p), 
with the Goldstone mode at the centre of the Brillouin 
zone, and a gapped branch (n_,p), with the lowest energy 
excitation at the X point [see Fig.j^c)]. The presence of 
the Goldstone mode indicates the breaking of a continu¬ 
ous SU(2) symmetry as expected for the correlated Ghern 
insulator ( |^ , see Sec. HA for details. Finally, it should 
be mentioning that the spectrum shown in Fig. [^is qual¬ 
itatively similar to the spin-wave excitation spectrum of 
the two-dimensional Mielkes model (flat-band limit) de¬ 
rived by Kusakabe and Aoki in the weak-coupling regime 
[see Fig. 1(a) from Ref. HH]. 


IV. FERROMAGNETISM IN A ELAT-BAND Za 
TOPOLOGICAL INSULATOR 

In this section, we study the flat-band ferromagnetic 
phase of a topological Hubbard model that preserves 
time-reversal symmetry. Due to the similarities with 
Sec. |HI[ here we just quote the main results and comment 
on the differences between flat-band Ghern and topolog¬ 
ical insulators. 


A. Time-reversal symmetric square lattice tt-Aux 
model 

Similar to Sec. Enl we consider N spinfull nonin¬ 
teracting electrons hopping on a bipartite square lattice 



K 


FIG. 4. (Color online) Dispersion relation (651 of t he ele¬ 
mentary excitations of the effective boson model (581 along 
paths in the Brillouin zone [Fig. [^c)] at the harmonic ap¬ 
proximation. Such a spectrum corresponds to the spin wave 
excitations of the flat-band ferromagnetic phase of the Ghern 
insulator (53l. 


where each sublattice, A and B, has Na = Nb = N 
sites. The Hamiltonian of the system T-Lq^^ is given by 
Eq. (0 but now we assume that 7(t) = —7(4-) = 1 for 
the nearest-neighbor hopping energy ([^. Such a choice 
implies that time-reversal symmetry is preserved, see Ap¬ 
pendix]^ for details. Bq^^ can be seen as two copies of 
the spinless tt-Aux modeP^ where electrons with spin f 
and j, are in the presence of opposite (fictitious) staggered 
±7r flux patterns, see Fig. Wa). 

In momentum space, assumes the form ([^ but 

with — h[_^, which is precisely a consequence of in¬ 
variance under time-reversal transformations. It is easy 
to see that Bq^^ can be diagonalized by the the Bogoli- 
ubov transformation 


•-kAt “ ^kdk^ + '^^kOjk^ 




(67) 




= u(lk4l +^-kcL. 




with the Bogoliubov coefficients Uk and Uk given by 
Eq. (BI I. Note that since 4 = 4ki the canonical trans¬ 


formation for spin f electrons differs from the one for spin 
j, electrons in contrast to the Ghern insulator discussed 
in Sec. |HI A where both transformations are equal [see 
After diagonalization, the Hamiltonian also 
(18), but now the 4x4 matrix 77'k 


Eq. ([^ 
assumes the form 
reads 


B'v = ( 


0 


-k 


( 68 ) 


with the 2x2 diagonal matrix h'^ giving by Eq. 
Similarly to the Ghern insulator discussed in Sec. I 


( 20 ) 
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the resulting band structure comprises two doubly degen¬ 
erate bands c and d, whose dispersion relations iOdd.aM 
are also given by Eq. (211 [see Fig. [^. As required by 
time-reversal symmetry® = ^^c/d-a-k- In par¬ 
ticular, for the TT-flux model , u}c/d,^,k. = ‘^c/d,iM- 

Again, the c and d bands are topologically nontrivial. 
Indeed, it follows from Eq. ( |2^ in combination with the 
fact that 7(t) = —7(4-) = 1 that the Chern numbers of 
the c and d bands are given 


Cf = -Cf = -1 


and 


q = -q = +1, 


i.e., as required by time-reversal symme¬ 

try. As a consequence, the charge Chern numberJi^ of 
the c and d bands vanish, 




:/d 


cf) 


= 0 , 


while the corresponding spin Chern numbers are nonzero. 


^cjd ^ 1 


c/d _ ^ 


Since c onserves the i-component of the total spin 

(see Sec. IIA), the Z 2 topol ogica l invariant for the c and 
d bands are simply given bjffEni 

t^c/d = mod 2 = 1, (69) 

implying that is a Z 2 topological insulator. Indeed, 

at half-filling (configuration not considered here), 71™^ 
should display a QSHE with the spin Hall conductivity 
= eCsl2Tr. 


B. Bosonization scheme for a flat-band topological 
insulator 


In this section, we introduce a bosonization scheme for 
a flat-band Z 2 topological insulator similar to the one for 
the fiat-band Chern insulator discussed in Sec. Im^ As 
shown below, the two bosonization schemes are quite sim¬ 
ilar, but there are important differences due to the fact 
that here time-reversal symmetry is preserved. Again, we 
focus on the nearly flat-band limit of the tight-binding 
model (^2 = at 1/4-filling {N^ = Nb = 

N). We restrict the Hilbert space to the lowest energy c 
bands and also assume that the ground state of is 


given by the ferromagnet state (24). 

Instead of the spin operator (251 at site i of the sub¬ 


lattice a, we now consider the following spin operator 
1 


^iab — d 




' ^ib c 


(70) 


with (a, 6) = {A,B) and {B,A). Using Eqs. (|^ and (26), 
the expression of the spin operators (70) in momentum 
space can be derived. In particular. 


^Cl,ab 'y A d— 

P 

P 


(71) 


wliere ± The spin operators 

projected into the c bands, i.e., the equivalent of Eq. (29), 
now read 


-^q.ab = -q)Cp-qTCp4, 

P 

-^q.ab = q)Cp-qqPt’ 


(72) 


where 


Cyis(p, q) — V—pBci^p and Gi 5 ^(p,q) — zi_p_|-qr?p. 

Again, we consider the following linear combination of 
the spin operators 

* 5 q.a = + (-l)“ 5 q,BA (73) 

with X = x^y^z and a = 0,1. In particular, we have 
^ ^ 9a ( P; q'j'^P 4 '’ 


'S'q.c = 

p 

where the gct(p,q) functions are now defined by 

5Q(P:q) “ p-|-qWp ( 1) "I*— p+q^p 


(74) 


(75) 


with a = 0,1 [compare Eqs. (33) and (75)]. 

Similar to the flat-band Chern insulators, the algebra 
of the projected spin and electron density operators is not 
closed. For instance, the equivalent of the commutator 
now reads 


['5'4;a>>§q'./3] [5a(-P + q'.-q)5/3(P,q')Cp-q-q'TCpT 


- 5a(-P:-q)'7/3(p- q,q')Cp_q_q/;Cp; .(76) 

The complete algebra of the projected spin and electron 
density operators can be found in Appendix [C| 

The construction of the boson operators ba in terms 
of the fermion operators c follows the same procedure 
outlined in Eqs. ([^-(43). Importantly, the spin opera- 

is 


tor that defines the particle-hole excitation [Eq. 
now given by Eq. (73). Again, we can define two sets of 


independent boson operators, namely 


bcB = ^5l5a(-P>q)< 

-Crv.o 


= — 


■p+qt'^qf’ 




(77) 


with a = 0,1, that obey the commutation relations (44). 


Here, the gQ(p,q) functions are given by Eq. (75) and 
T^a.q = ^g*c.i-P + ^, q)5a (p, q) 


(78) 




















11 


[see Eq (Cl) for the expression of the function in 
terms of the coefficients It is worth mentioning that 

although for both Chern and Z 2 topological insulators the 
bosons ba are linear combinations of particle-hole pair 
excitations, for the former the particle and the hole are 
on the same sublattice [see Eqs. ( |M| ) and (43)] while, 
for the latter, the particle and the hole are on different 
sublattices [see Eqs. (73) and (ff^j.^ 


Einally, the electron density operator (36) projected 
into the c bands [the equivalent of Eq. (37 ijiiow reads 


Pa a (k) ^ ^ Fj (t(Pj 

P 


p-ko-'^pt 


(79) 


with 


G^^(p,k) = V 


p-k^pi 


GBt(P’k) = u 


p—k^P ’ 


(80) 


Gai{p, k) = -u.p+kP* p, Gb lip, k) = u-p+kudp. 


Eollowing the procedure outlined in Eqs. (46|-(52), the 


bosonic representation of (79) can be easily derived. It is 


also given by Eq. (52) but with the Gapaa{x,y) function 
now given by Eq. (C5). Here, the equivalent of Eq. (51) 
reads 


p, q) = 7 ^— 9a{-p + q, q)- 


F, 


(81) 


Q;,q 


C. Topological Hubbard model II 


In this section, we consider a correlated topological in¬ 
sulator on a square lattice described by the Hamiltonian 

T-Lz2 =T~Lq^^ FT-Lu, (82) 

where is the square lattice tt-Hux model discussed 

in the nearly flat-band limit (^2 = 


Sec. 


IV A 


and 77(7 is the repulsive on-site Hubbard term ( |54| ). 
Again, 1/4-filling is assumed. The topological Hubbard 
model ( [^ has rece ntly b een discussed by Neupert et 
Similarly to Sec HI B we now apply the bosoniza- 
tion formalism introduced in Sec. jlVBj to study the flat- 
band ferromagnetic phase of the Hamiltonian (821. 

Eollowing the lines of Sec. [m^ the first step is to 
project the Hamiltonian (82) into the lowest-energy c 
bands, i.e. 


7722 —>■ 7722- 


We then map the projected Hubbard model 77 22 to an 
effective interacting boson model 77^^. We find that 77^^ 
has the same form as the boson Hamiltonian (58) but 
with the Gapaaix^y) function given by Eq. (C5). 


Within the harmonic approximation, the effective bo¬ 
son model 77f^ can be diagonalized and it assumes the 
form (64). The dispersion relation of the bosons a± is 


equal to (65), i.e. 


nil = ±Ap 


f2 _ ^10^:01 

p p p 


(83) 



K 


FIG. 5. (Color online) Dispersion relation (831 of the elemen¬ 
tary excitations of the effective boson model 77f ^ along paths 
in the Brillouin zone [Fig. [^c)] at the harmonic approxima¬ 
tion. Such a spectrum corresponds to the spin wave excita¬ 
tions of the flat-band ferromagnetic phase of the Z 2 topolog¬ 
ical insulator (821. 


but with the Gapanix^y) function given by Eq. (C5). 


Simil arly t o the correlated Chern insulator discussed 
in Sec. HIC the ground state of 77is the vacuum 
state for the bosons a±, i.e., the ferromagnetic state |EM) 
[Eq. (|24])b and the excitation spectrum [the dispersion 
relation (83) of the bosons a±] of such a flat-band fer- 
romagnet has two branches, see Eig. However, in 
contrast to the spin-wave spectrum of a Chern insula¬ 
tor (Fig. 1^, here both branches are gapped at zero wave 
vector. Differently from the flat-band Chern insulator 
(53), where a continuous SU(2) symmetry is broken, the 


ground state of the correlated topological insulator pre¬ 
serves the U(l) spin rotation symmetry of the Hamil¬ 
tonian (82) (see Sec. HA for more details). As a con¬ 
sequence, a Goldstone mode is absent in the excitation 
spectrum. Again, the energy scale of the spin-wave exci¬ 
tations is given by the on-site repulsion energy U. The 
excitation gap A = = 0.3134 7/, which is at the 

centre of the Brillouin zone, agrees with the exact diag- 
onalization data from Neupert et al, A « 0.30 7/.!^ The 
fact that the spin-wave excitations remain gapped at all 
wave vectors corroborates that the ground state is indeed 
given by our reference state (24) and not by an in-plane 


(W-type) ferromagnetic state. 


Interestingly, since a phase with ferromagnetic long- 
range order sets in, time-reversal symmetry is sponta¬ 
neously broken. As shown in Ref. [THl the ferromagnet 
ground state also displays an IQHE with the Hall con¬ 
ductivity axy = jh. 
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V. DISCUSSION 


Although we have focused our discussion on the square 
lattice TT-flux model, the bosonization formalisms for flat- 
band Chern and topological insulators, respectively in¬ 
troduced in Secs. Ell and |IVB[ are rather general. In 
principle, they can be employed to study the flat-band 
ferromagnetic phase of a topological Hubbard model 
whose single-particle term assumes the 4x4 matrix form 
Q, such as the Kane-Mele-Hubbard model without the 
Rashba spin-orbit coupling.l^In this case, spin is a good 
quantum number. Once the coefRcients i?i.k [Eq. 0] 
of the model are identified, the corresponding effective 
(interacting) boson model is easily determined since the 
coefficients (61) and (62) of the boson model are written 


in terms of the = Hi_k/|Bk| functions (see Appen¬ 
dices E] and [^ . One important point is to verify whether 
the condition ( [4^ holds, i.e., if it is possible to define 
two sets of independent boson operators. It would be 
interesting to see whether the bosonization scheme could 
be extended to the following cases: (i) Four band mod¬ 
els where spin is not conserved. It would allow us to 
consider, e. g., the Kane-Mele-Hubbard model with a 
Rashba coupling.!^ (ii) Six band models where the single¬ 
particle term assumes the form 0 but with being a 
3x3 matrix. Two examples are the tight-binding model 
on the kagome lattic^ and the three-orbital square lat¬ 
tice model.^ 

One interesting feature of the bosonization scheme de¬ 
veloped here is that it allows us to analytically deter¬ 
mine the spin-wave excitation spectrum of a flat-band 
ferromagnet with topologically nontrivial single-particle 
bands. At the moment, it is not clear how to compare the 
results derived in Secs. |HI C| and |rV^ with other approx¬ 
imation schemes. For the 2DEG at filling factor v = 
the noninteracting term of the effective boson model de¬ 
rived within the bosonization formalisrrfi^ describes the 
magnetic exciton excitations of the quantum Hall ferro¬ 
magnetic ground state. The energy of the bosons are 
exactly equal to the magnetic exciton dispersion rela¬ 
tion derived by Kallin and HalperiiP^ within diagram¬ 
matic calculations. Such a diagrammatic formalism is in¬ 
deed equivalent to the so-called time-dependent Hartree- 
Fock approximation.!^ Although for flat-band Chern and 
topological insulators only an approximated bosonic rep¬ 
resentation for the electron density operator can be de¬ 
rived [Eq. ([^], we expected that, due to the analogy 
with the 2DEG at = 1 (see discussion at the end of 
Sec. HI B), the spin-wave dispersion relations determined 
in Secs. HI C| and |IV C agree with a time-dependent 


Hartree-Fock analysis of the topological Hubbard models 
(pi and (pi. 


A second interesting aspect of the bosonization scheme 
for flat-band Chern and topological insulators is that it 
provides an interaction between the bosons (spin-waves). 
A detailed study of the consequences of the spin-wave- 
spin-wave coupling is beyond the scope of the present 
paper. It would be interesting, for instance, to verify 


whether the boson-boson interaction (60) yields two-spin- 
wave bound states and whether such bound states are re¬ 
lated to possible topological excitations of the flat-band 
ferromagnetic state (24). Such an analysis is motivated 
by the fact that for the 2DEG at filling factor v = 
the boson-boson coupling derived within the bosoniza¬ 
tion formalisnJi^ gives rise to two-boson bound states 
that are related to skyrmion-antiskyrmion pair excita¬ 
tions (the charged excitation of the 2DEG at = I 
is described as a topological excitation, quantum Hall 
skyrmion, of the quantum Hall ferromagnetic ground 
state). This set of results allows us to properly treat the 
skyrmion as an electron bound to a certain number of bo¬ 
son (spin-waves) We defer the analysis of the effects of 

the boson-boson interacting (60) to a future publication. 


VI. SUMMARY 

In this paper, we have considered the flat-band fer¬ 
romagnetic phase of a correlated Chern insulator and a 
correlated Z 2 topological insulator and analytically cal¬ 
culated the corresponding spin-wave excitation spectra. 
In particular, we have considered two variants of a topo¬ 
logical Hubbard model, namely, one with broken time- 
reversal symmetry (Chern insulator) and another one 
invariant under time-reversal symmetry (topological in¬ 
sulator). In both cases, the single-particle term is the 
square lattice tt-Hux model with topologically nontrivial 
and nearly flat bands. The spin-wave dispersion rela¬ 
tion has been determined within a bosonization scheme 
similar to the formalisirfi^ previously developed for the 
2DEG at filling factor v = 1. Here, we showed that 
the formalisrrPS can indeed be generalized for flat-band 
Chern and topological insulators. 

The investigation of the spin-wave excitation spec¬ 
trum indicates the stability of the flat-band ferromag¬ 
netic phase. Generically, we obtain two spin-wave ex¬ 
citation branches due to the two-atom basis of the lat¬ 
tice model. We find that the correlated flat-band Chern 
insulator has a gapless excitation spectrum: the Gold- 
stone mode is associated with a spontaneous SU(2) spin- 
symmetry breaking. For the correlated flat-band topo¬ 
logical insulator with preserved time-reversal symmetry, 
the excitation spectrum is gapped since the flat-band fer¬ 
romagnetic ground state preserves the U(l) spin rota¬ 
tion symmetry of the Hamiltonian. Moreover, within the 
bosonization scheme, we find a spin-wave-spin-wave in¬ 
teraction. Due to the analogies with the 2DEG at filling 
factor = 1, we expect that such coupling may give rise 
to two-spin-wave bound states. 
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Appendix A: Symmetries of the yr-flux model: 
time-reversal 


Tr^^T ^ and Tr^^T ^ we see that only 

the Hamiltonian (A3) is invariant under time-reversal. 


In this section, we discuss the behavior of the nonin¬ 
teracting fermion model ([^ under time-reversal. 

The time-reversal operator T is defined as 

T = i{<Ty ® I) K, (Al) 

where K denotes complex conjugation and I is the 2 x 
2 identity matrix. Invariance under time reversal, i.e., 
[HoyT] = 0 , implies tha 1 p^ 

T^kT-i = H-u, 

where "Hk is the matrix ([^. Since 


Appendix B: Details of the bosonization scheme for 
flat-band Chern insulators 

In this section, we quote the expressions in terms of 
the coefficients Bi y^ of some quantities related to the 
bosonization scheme discussed in Sec. nni 

Let us first consider the diagonalization of the nonin¬ 
teracting Hamiltonian 0 - It is useful to write the coef¬ 
ficients Uk and Uk of the Bogoliubov transformation 0 
as 


THkT-^ 


hi* 0 
0 hi* 


invariance under time-reversal implies that hi = as 
mentioned in Sec. El 

Alternatively, we can follow Lu and Kan^^ and write 
the matrix "Hk in terms of the five 4x4 Dirac matrices 


pl,2,3,4,5 ^ ^ ^ Ty,ax 0 Tz,ay (g) g) Tz) , 

and their ten commutators T*-! = [T®,r^]/(2i). Here, 
^x,y,z and Ty j, 2 are 2x2 Pauli matrices respectively re¬ 
lated to spin and sublattice. The Dirac matrices obey 
the Clifford algebra r^T-l -|- T^T-^ = 26ijl. For the Chern 
insulator discussed in Secs. |H] and 

Hk = SpkFi + Ba.kF^ + Ha.kTi" (A2) 


HI A 


while, for the topological insulator considered in Secs. |TT] 
and IIV Al 


Uk = exp(-|-z(/)k/ 2 ) cos(dk/ 2 ), 

Wk = exp(-z(/)k/2) sin(6»k/2). (HI) 

Due to the form of the matrix hi [see Eq. @] , it is in¬ 
teresting to introduce the following relations between the 
(f>y_ and 0k functions and the coefhcients k [Eq. 0 ]: 

Bi.k = sin 0 k cos 0 k, ^ 2 ,k = sin 0 k sin 0 k, 

Ba.k = cos 0k, 

sin 0 k = J Bfy. + i?| 1 ^, tan 0 k = (B 2 ) 

^ ’ ’ 5i,k 

where Bk = Bk/|Bk|. It then follows that 

l^kP = 2 , kkP = 2 , 

“k^k = ^ (A,k + *. 62 ,k) ■ (B3) 


Hk = BpkFi + + Ba.kF^^ (A3) 

Since k _k, 

^pip -—1 _ f -fF , i 1 , 

" 1 -r®, * = 2,3,4,5, 


With the aid of Eqs. (B1)-(B3), it is easy to show that 
the Hamiltonian (18) is the diagonal form of ([^ . 

Concerning the algebra of the projected spin and elec¬ 
tron density operators, in addition to the commutator 
(Ml, it follows from Eqs. (IM) and (IM) that 


[pac7(k), S'+„] = Yi [^<T,tG'a(p - q, k)gQ,(p, q) - i5a,;Ga(p, k)ga(p - k, q)] (^4) 

P 

[parT(k), ^ [^<T,;Ga(p - q, k)g„(p, q) - Sa^^Gaip, k)ga(p - k, q)] c^_q_k;Cpt, (B5) 


[parT(k), Pba' (q)] = ^^,<7' Y k)Gb(p, q) - Ga{p, k)Gb(p - k, q)] cj 


p-q-k(T®'P<^’ 


(B 6 ) 


with a,b = A, B, a = 0,1, and the Ga(p, q) and ga{p, q) functions respectively given by Eqs. (30) and (33). 

Within the approximation (40), the commutator (38) assumes the form (42). With the aid of Eq. (B3), it is possible 
to show that {a ^ (3) 


Y ‘l)00(P - q: -q) = ^ XI (^3.p + B3,p_qj -f *(-!)“ (i3l,p-q.B2,p - H 2 ,p-qBl,p^ 


(B7) 
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Since for the tt-Aux model, [see Eq. 0] it is then easy to show that Eq. ( |B7[ ) vanishes. 

Similar to Eq. (B7), it is possible to write the function in terms of the coefficients Using Eq. (B3), 
show that 


we 


~ 2 -®3,pE3^p-q( 1) ^.Bl^pi?l^p_q + B2^pB2,p-cij 


(B8) 


The (approximated) boson representation of the projected electron-density operator is given by Eq. (52) where the 
Qapaa (k, q) function is defined as 


0a/3at(k,q) = ^-VGa(p,k)g„(p-k,q)g/ 3 (p-k-q,-k-q), 

^a,q^/3,k+q p 

0a/3a;(k, q) = +-—^^ Ga{p - q, k)g„(p, q)ff/3(p - k - q, -k - q). 
^a,q^/9,k+q 


(B9) 


Again, using Eqs. (B3|, we find after some algebra that 

1 


1 


0a /3 a T (k, q) = - q [<5a,A + (- 1)“+'^] 


-^a,q-^/3,k+q 


Ei-3(-ir^3.p 


TE3 p_q.B 3 p_|_k T B^ p—f^B^ p -)- .B3 p_|_kE 3 ^p ( 1) B3 p_q.B 3 P.B3 p-i-k 

—l)“[i?l,p-q.Bl,p + E2,p-qi?2,p + * (” 1)“ (El p_qi?2,p — E2,p-qEl,p)] [1 — (—l)“E3,p+k] 

+ (~l)^[El,p-qEi p+k + E2^p_qi?2,p+k ~ *(~ (El^p_qE2.p+k ~ E2,p_qEl ^p+k)] [1 ~ (~1)“E3 p] 

+ (-l)“+^[Ei,p+kEi,p + 4^p+kE2.p -b z(-l)“(Ei,p+kE2^p - 4.p+kEi,p)][l + (-l)“4,p-q]. (BIO) 


The expression of ^^^^^(k, q) can be derived from Eq. (BIO) using the fact that 0a/3a^(k, q) = —k, —q). 

Finally, the bosonic representation of the projected single-particle Hamiltonian (57). The first step is the calculation 
of the commutator 

[Eo, foLq] = ^"i^’^^ Cp-qtCpt = (‘^c.p-q ” Wc,p) g/3 (k, P, q)6^^k- (BH) 


',k,p 


Here, the second equality follows from Eq. (47). Note that the expansion 

5a(p,q^ 


T~('Q,B — ^ ^ ^ ^ (^c,p—q ^c,p) 

q:,/ 3 k.p,q 




E/3(k,p,q)6^_i^6„.i 


(B12) 


of "Ho in terms of the boson operators ba satisfies the commutator (Bll). Keeping only the term k = q in the 
momentum sum above and using Eq. (51), we have 


Eo.s = Eo-b ^ ^a;a/3(q)5^_q5, 


Q!,q5 


(B13) 


a.p q 


where Eq = 2.44A^ is a constant related to the action of T-Lq in the reference state (24) and 

1 


Wa/3(q) = 


Ea,qE^ q 


^ [ (a^e,p—q Wc,p) ffa(P) q)<?/3 (p qj q)■ 


(B14) 


In the flat-band limit, a;Q,^(q) = 0 since Wc^p = 0. In the nearly flat-band limit considered in Sec. HIC the coefficient 
Wa/3(q) can be nonzero. However, for the tt-Aux model, it is possible to show that Wa/3(q) vanishes due to the fact 
that Ej k “ E^ _k. 
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Appendix C: Details of the bosonization scheme for flat-band topological insnlators 


Similar to the previous section, we here provide the expansion in terms of the coefficients Bi \^ of some functions 
related to the bosonization formalism introduced in Sec. lYl 


The function (78) reads 

^ E [(-!)“(! - 4,pB3.-p + q) + Sl,pBl,_p+q + 4,p 


F. 


2,-p+q 


(Cl) 


The algebra of the projected spin and electron-density operators. From Eqs. (74) and (791, we show that, in addition 
to the commutator (76), the following commutation relations hold 


[Paa(k), S'+„] = ^ [5cr,^Gaa{v “ Q, (-p, -q) - 5c,,iGaa{v, + k, -q)] (C2) 

P 

[Paa(k),5q^] = ^ ^ (p - q, k)^^ (p, q) - 5,T.tG'aa(p,k)g„(p- k,q)]cj,_q_i^^Cp^, (C3) 

p 

[pa a (k) , Pb a' (q)] = E (P “ a' (P, q) - Ga ,7 (P, k)Gb < 7 ' (p “ k, q)] ^Cp (C4) 


with a,b = A, B, a = 0,1, and the gaip, q) and Gaa{p, q) functions respectively given by Eqs. (75) and (80). 


The Gafiaa-O^, q) function, which is related to the bosonic representation of the projected electron-density operator 

1 


(79), is defined as 


Ga P a t(k, q) — 


^a/3a^(k,q) — -)- 




1 


.^a,qF,3,k+q 


E k)5a(p - k, q)g^(-p -b k -|- q, k + q). 


E k)g„(p, q)g^(-p -b k -b q, k -b q) 


(C5) 


and, in terms of the coefficients Bi q, it reads 


G.patiK q) = - + Sa,Bi-ir+^]-^ - E(-l)“[l - (-l)“^3.p 


Fa,qFp,k+q p 

-^3, —p+q'^3,p+k .^3, —p-|-q-^3,p T -^3,p+k.^3,p ( 1) -^3, —p-fq-^3,p-^3,p+k] 

+ [-Bl,_p+qi3l,p -b F2,-p+qF2,p -b *(— 1)“(Fl,_p+qi?2.p ” F 2 ,_p+q.Bl^p)] [1 — (—l)“i33_p+k] 

+ (~l)'^~''^[-®l,-p+q-®l,P+k + F2,_p+qi32,p+k — *(—l)“(.Bl,-p+q.B2^p+k ” .B2,_p+q.Bl,p+k)] [1 ” (” l)“.B3^p] 

+ ( — p+kFi^p -b .B2,p+kF2,p + *( —l)'^(.Bl,p+k.B2,p — ^2,p+kFi^p)][l — ( —l)“.B3__p+q] (C6) 

The expansion for q) follows from Eq. (C6) using the relation 5Q,^a^(k, q) = —(—l)“+^5a^at(k, q). 
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